Geometrically induced modification of surface plasmons 
in the optical and telecom regimes 
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We demonstrate that the introduction of a subwavelength periodic modulation into a metallic 
structure strongly modifies the guiding characteristics of the surface plasmon modes supported by 
the system. Moreover, it is also shown how a new type of a tightly confined surface plasmon polariton 
mode can be created by just milling a periodic corrugation into a metallic ridge placed on top of a 
metal surface. 



Photonics based on the exciting capabilities of surface 
plasmon polaritons (SPPs) has become a very active area 
of research during the last decade [TH3]. Strong local- 
ization of electromagnetic (EM) fields and the building 
up of ultra-small SPP-based waveguides [4j [5] have been 
achieved thanks to the subwavelength (the term "sub- 
wavelength" refers to the vacuum wavelength) nature of 
the SPP fields, thus enabling a great variety of applica- 
tions in optics [6j [7] . 

However, in order to fulfill the potentialities of light 
guiding based on SPP excitation, it is convenient to 
search for effective ways to control and tune the prop- 
agation characteristics of SPP modes. The concept 
of geometrically-induced SPPs [5HT3] (also named spoof 
SPPs) has proven to be very powerful in tailoring the 
dispersion relation of the propagating surface EM modes 
in the microwave and terahertz regimes. The aim of this 
Letter is to transfer the spoof SPP concept to the opti- 
cal and telecom ranges. We will demonstrate how the 
dispersion relation of propagating SPP modes can be 
tuned by introducing a periodic modulation in the metal 
surface in a length scale much smaller than the wave- 
length (A). Even more, we will show how this modula- 
tion could also "build-up" tightly confined SPP modes in 
structures where these modes were not supported in the 
non-corrugated case. 

First we consider a SPP mode that is bound to and 
propagates along a finite V-shaped groove milled in 
a metal film, the so-called channel plasmon polariton 
(CPP). The propagation characteristics of these CPP 
modes have been studied both theoretically [T4HT6] and 
experimentally [17] during the last years, and even pro- 
totype CPP-based photonic circuits have been already 
fabricated [T8] . In Fig. [TJa) we plot the dispersion re- 
lation (frequency versus parallel momentum) of the fun- 



damental CPP mode for a V-groove whose geometrical 
parameters are taken from experiments [19]. The metal 
considered is gold and its frequency-dependent dielectric 
function is taken from [20]. All the numerical calcula- 
tions presented in this paper have been performed using 
a Finite Element Method (FEM). Due to the finite depth 
of the grooves, this CPP mode presents a cutoff wave- 
length that, for the chosen set of geometrical parameters 
(see caption of Fig. [I]), appears at around 1.3 /im. Grey 
curves in Fig. [TJa) illustrate the modification of the CPP 
dispersion relation induced by a sub-A periodic modula- 
tion in the V-groove. The main effect of the corrugation 
is to shift the cutoff wavelength of the CPP mode to 
longer wavelengths, reaching a value of 1.6 /im for the 
shortest period analyzed, d = 25 nm. Also, as its disper- 
sion relation departs more from the light line, the CPP 
for the corrugated V-groove becomes more localized than 
the pure one. 

Our finding that cutoff wavelength shifts to longer 
wavelengths in corrugated V-grooves could explain why 
in the experiments there exists a propagating CPP mode 
in the wavelength range between 1424-1640 nm [17], de- 
spite the fact that calculations for a non-corrugated V- 
groove with the same geometrical parameters (6 = 25 
degrees) predict a cutoff wavelength of 1440 nm [TB] . 
Scanning electron microscope images reveal the presence 
of a weak periodic modulation in V-grooves fabricated 
with focused ion beam techniques. When a shallow cor- 
rugation with a period d = 25 nm is now introduced in 
the V-groove, our calculations show that the cutoff wave- 
length moves from 1440 nm to 1750 nm, larger than the 
wavelength range analyzed in the experiments. 

A much stronger effect associated with a longitudinal 
sub-A periodic corrugation is seen for another type of 
SPP modes: a slot SPP mode that propagates along a 
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FIG. 1. (a) Dispersion curves for CPPs. Red and blue curves 
represent those for a CPP mode of the V-groove without cor- 
rugation for two different angles, = 30° and = 15°. The 
depth of the V-groove is h = 1.1 /im and the radii of curvature 
at the edges are r = 10 nm and R = 100 nm. Grey lines dis- 
play the dispersion curves for CPP modes of the corrugated 
V-groove (0 = 30°) where the depth of the grating is fixed at 
30 nm and three different periods are studied: d = 400 nm, 
d = 100 nm and d = 25 nm. Inset shows the geometry, (b) 
Dispersion curves for slot waveguide modes. Red and grey 
curves are those for the slot without grating. Blue curve is 
the dispersion curve for a corrugated slot. The period of the 
grating is d = 100 nm. 

gap carved between two metal plates. Here we consider a 
waveguide structure that is composed of a slot perforated 
on a thin silver film deposited on top of a glass substrate 
(e = 2.25), see inset of Fig.[TJb). The main geometrical 
parameter that controls the dispersion relation of these 
modes is the width of the slot, w. Figure [TJb) shows 
two different non-corrugated structures, w\ = 75 nm (red 
line) and = 135 nm (grey line). The geometry of the 
corrugation is chosen so the distance between two deepest 
points, W, is equal to W2 whereas the minimal distance 
is equal to w\, see inset of Fig. [TJb). The introduction of 
a periodic modulation has a dramatic effect on the dis- 
persion relation. One could naively expect that the re- 
sulting dispersion curve would be located between those 
for the non-corrugated cases. However, the band for the 
corrugated slot departs strongly from those two curves, 



resulting in a much longer cutoff wavelength and larger 
confinement. Then, our results clearly show that a sub- A 
periodic modulation could also be used to improve the 
guiding properties of slot waveguide modes. 

Now we demonstrate that a sub- A periodic modulation 
could indeed create a SPP mode in structures where, 
without corrugation, laterally confined SPP modes are 
not supported. As an example we consider a metallic 
ridge (height h and width L) placed on top of a substrate 
made of the same metal, see Fig. |2|b). This structure 
does not support the propagation of SPP modes transver- 
sally localized. In Fig. [2|a) we render the dispersion 
relations of the geometrically-induced SPP modes that 
emerge when a sub-A periodic modulation is introduced 
into a gold ridge. The six curves correspond to different 
values of the modulation depth, g. As clearly seen in the 
figure, even the weakest modulation (g — 20 nm) is able 
to create a laterally confined SPP mode (i.e., the disper- 
sion curve is lower than the SPP-curve for the flat metal 
surface). When the depth is enlarged, the dispersion 
relation further departs from the light line, increasing 
the mode localization. Accompanying this movement, 
the cutoff frequency shifts to lower frequencies. The in- 
crease in the mode localization also affects the propaga- 
tion length of these geometrically-induced SPP modes. 
For a fixed wavelength, the transversally confined mode 
decays faster for a larger g, see Fig. [2|c). 

In the case where the grating depth is equal to the 
height of the slab, the geometry resembles a ID chain of 
metallic box-shaped particles placed on top of a metal 
film. From now on, we name the mode supported by this 
structure as Domino Plasmon-Polariton (DPP). This 
DPP mode has a high electric field localization near the 
top part of the domino structure, see Fig. [3] In this 
figure the horizontal slice renders the electric field inten- 
sity evaluated in a xz-plane that is parallel to the metal 
substrate and located 5 nm above the domino's top face. 
The intensity also presents a strong subwavelength con- 
finement in the transversal direction. Regarding the vec- 
torial nature of the EM-fields associated with a DPP, the 
electric field has mainly y and z components (see yellow 
lines in the vertical plane of Fig. [3| whereas the magnetic 
field has predominant x and z components (see blue lines 
in the horizontal plane of Fig. [3|. 

It is worth discussing the differences between the DPP 
modes described above and the plasmon modes sup- 
ported by ID arrays of metal nanoparticles placed on 
top of a dielectric film [2TJ [22]. In this last case, the 
near-field coupling between the localized plasmon modes 
of the metal nanoparticles leads to the formation of a 
very flat plasmon band characterized by a deep sub-A 
confinement but short propagation length. However, in 
the case of DPP modes, the presence of the metal sub- 
strate results in the emergence of a "polaritonic" part in 
the dispersion relation that runs close to the SPP band of 
the flat surface [see Fig. |2ja)]. Accordingly, DPP modes 
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FIG. 2. Creation of a geometrically-induced ridge SPP mode, (a) Dispersion curves for SPP modes running on a corrugated 
ridge with different modulation depths, (b) Geometry of the ridge structure with a modulation: corrugation period d = 75 nm, 
width of the ridge L = 37.5 nm, height of the ridge h = 120 nm and the grooves width a = d/2 = 37.5 nm. (c) Corresponding 
propagation lengths for the cases analyzed in panel (a). 




FIG. 3. Electric and magnetic fields associated with a domino 
plasmon polariton mode. Yellow (blue) lines represent the 
electric (magnetic) vector field. Electric field intensity map is 
evaluated on a horizontal plane located on top of the domino 
structure. The geometrical parameters of the domino struc- 
ture are: modulation period d = 75nm, lateral width of the 
ridge L = 37.5 nm, height h = 120 nm and the grooves' width 
a = d/2. The operating wavelength is A = 1.5/zm. The defi- 
nition of the cartesian axes is also depicted. 

operating within this polaritonic regime posses a very 
long propagation length but, as expected, are much less 
confined than the modes supported by a chain of metal 
nanoparticles. 

The propagation properties of a DPP mode depend on 
the lateral size (L) of the structure. By varying this di- 
mension we could find optimal properties of the DPP 
mode, like enhanced field confinement and/or longer 
propagation length. Figure [4^a) shows the DPP disper- 



sion curves for four different values of L. As this width 
is reduced from L = oo to L = O.bd = 37.5 nm, the 
cutoff frequency increases and, for a fixed frequency, the 
wavevector is smaller for narrower dominos. Accordingly, 
the DPP propagation length for very narrow dominos is 
larger than that for a wider one (see Fig. Qb)). Interest- 
ingly, even for a domino width of only L = d/2 = 37.5 nm 
the propagation length at A = 1.5/im is around 80 /im, 
a value that is larger than those predicted for other 
types of SPP modes that present sub- A confinement like 
CPPs [16] or wedge plasmon polaritons [23]. When the 
metal behaves as a perfect electrical conductor, our cal- 
culations show that the dispersion relation of the DPP 
mode is almost independent of L, paving the way to sev- 
eral interesting applications for the use of DPPs in the 
terahertz regime [24] . 

Figures [4^ c-e) render the electric field intensity distri- 
bution for three DPPs in a plane perpendicular to the 
propagation direction located at the center of the gap 
between two dominos. The lateral sizes analyzed are 
L = 0.5d, L = 3.5d and L = lOd and the operating 
wavelength is A = 1.5/im. The color scale is the same 
for the three distributions. As expected from the dis- 
persion curves, the DPP mode is more confined for the 
wider domino, L = lOd. When the domino width is 
decreased, the EM energy is carried within a bigger vol- 
ume. This is illustrated in Figs. [4^c-e) by representing 
the energy isocurves for the three different dominos. The 
percentage associated with each isocurve measures the 
amount of EM energy (with respect to the total one) car- 
ried out through the area inscribed into the curve. This 
is the usual way to define the modal size of a waveguide 
mode [25]. However, a closer look at the intensity distri- 
butions shows that the electric field in the case L = 3.5d 
presents a stronger field localization than the L = lOd 
or L = O.bd cases. Nevertheless, a comparison of electric 
field distributions with the operating wavelength (white 
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FIG. 4. Dependence of the DPP propagation characteristics with the domino width, (a) Dispersion curves of the DPP modes 
for four different widths L of the structure. The geometrical parameters of the domino structure are the same as in previous 
figures, (b) Propagation length versus free-space wavelength for the four structures analyzed in panel (a), (c-e) Electric field 
intensity distributions evaluated at A = 1.5 /xm in a plane perpendicular to the propagation direction for three different values 
of L. The blue lines show several energy isocurves: the percentage measures the amount of energy (with respect to total one) 
carried through the area bounded by the corresponding isocurve. 



bars) in Fig. 4 illustrates the strong sub-A transver- 
sal confinement associated with the propagation of DPP 
modes. 

In conclusion, we have shown that the application of 
the spoof plasmon concept in the optical and telecom 
regimes allows tailoring of the guiding properties of SPP 
modes. Moreover, we have also demonstrated that the 
same approach leads to the emergence of guided modes in 
geometries where conventional SPPs are not supported. 
We have illustrated this finding by analyzing the so-called 
Domino Plasmon Polaritons, which exhibit an excellent 
trade-off between lateral confinement and propagation 
length. 
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